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摘 要 : 基于 图 形 学 分 析 技术 ,采用 追踪 算法 .矢量 分 析 法 ,利用 常规 MICAPS 数据 、NECP 逐日 再 分 析 资 料 和 EC 


细 网 格 资料 ,对 我 国 西北 地 区 致 时 天气 系 统 进行 识别 试验 。 
踪 算 法 能 够 成 功 识别 高 空 槽 ,高空 脊 等 主要 系统 ,并 可 得 到 其 曲率 特征 ;@) 利用 该 方法 识别 天 气 尺 度 系 统 , 系 统 中 


通过 试验 ,得 到 如 下 结论:Q 依据 天 气 学 定义 ,利用 追 


心 识 别 率 高 达 90% 以 上 ,平均 误差 小 于 0.5" ,而 高 空 模 识 别 准 确 率 为 74% ~87% ,平均 识别 误差 为 0.7° ~1.4°。 
利用 风 场 进行 订正 后 ,准确 度 明显 提高 误差 有 所 减 小 ,基本 可 以 满足 干旱 预测 、 预 报 业 务 需求 。 


关键 词 : 天 气 系统 ; 图 形 化 ; 自动 识别 ; 追踪 算法 


近年 来 全 球 气候 变 暖 日 益 严 重 ,全 球 干 旱 问 题 
日 益 罕 出 。 干 旱 是 西北 地 区 自然 环境 的 最 大 特 
点 ,其 分 布 着 全 国 85% 的 干旱 或 半 干 旱 土地 面 
只 。 人 研究 导致 西北 干旱 气候 的 主要 天 气 系统 ,分 析 
其 本 身 特 征 ,预测 其 变化 趋势 ,对 西北 地 区 的 未 来 发 
展 具 有 十 分 重要 的 现实 意义 。 目 前 ,在 天 气 预报 业 
务 系统 中 ,通过 影响 天 气 系统 的 类 型 判 识 . 位 置 判 
识 .强度 判 识 ,根据 概念 模型 进行 预报 推理 ,是 最 为 
基本 的 预报 方法 之 一 。 影 响 天 气 系统 识别 的 基本 业 
务 需求 为 :确定 模 、 峭 及 其 位 置 ,闭合 高 , 低 系统 中 心 
位 置 与 强度 ,对 指定 地 点 的 相对 距离 与 前 后 位 置 , 系 
统 的 移动 路 径 及 发 展 趋势 等 。 

常规 的 天 气 图 分 析 重 视 的 是 高 度 场 和 风 场 ,但 
已 不 能 满足 日 益 精细 化 的 预报 需求 。《 现 代 气 象 预 
报 业 务 发 展 规划 (2016 一 2020 年 )》 中 提出 要 发 展 
“多 源 数 据 融 合 与 再 分 析 技 术 ”“ 重 点 组 织 高 影响 
天 气 监测 预警 预报 技术 等 关键 领域 研发 。 

精细 天 气 图 分 析 要 求 一 个 高 度 组 织 化 的 分 析 方 
法 (MICAPS 综合 图 分 析 ) ,根据 不 同情 况 将 分 析 数 
据 和 预报 参数 放 在 同一 张 图 上 。 预 报 员 根据 综合 图 
显示 的 不 同情 况 , 利 用 主观 图 形 识别 能 力 进行 识别 ， 
判断 不 同情 况 可 能 造成 的 不 同 天 气 现象 。 因 此 ,能 
否 基于 图 形 学 方法 来 实现 天 气 系统 的 自动 识别 , 减 
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少 预 报 员 手 工分 析 天 气 系统 的 工作 量 , 提高 精细 化 
预报 能 力 就 成 为 近 几 年 新 的 研究 课题 -2 。 

近年 来 , 随 着 电子 计算 机 性 能 的 提高 ,计算 机 图 
形 识别 技术 也 取得 飞速 发 展 ,而 其 成 果 多 应 用 在 工 
商 和 医学 “~ ,如 条 形 码 的 识别 53) .汉字 笔画 自动 
识别 "建筑 图 纸 识别 “路 、 人 脸 识 别 系统 "等 
方面 。 在 地 球 科学 研究 方面 ,图 形 识别 技术 也 在 逐 
渐 发 展 ,如 利用 计算 机 图 形 识别 技术 与 CIS 技术 相 
结合 ,进行 地 理 信息 的 绘制 与 识别 ;利用 图 形 
识别 技术 进行 沙 洁 天气 监测 识别 -3 或 卫星 云图 、 
多 普 勒 雷达 图 的 系统 识别 5 -2 。 近 几 年 ,也 有 一 些 
学 者 开展 了 天 气 系统 自动 识别 技术 的 研究 ,但 大 部 
分 都 是 从 数学 角度 进行 计算 识别 下-"?。 如 吉美 英 
等 5 利用 读 取 欧 洲 高 度 场 格 点 数据 ,运用 统计 识别 
法 和 模糊 识别 技术 ,对 高 低压 中 心 和 覃 线 进 行 分 析 
识别 。 胡 文 东 等 “中 读 取 MICAPS 格 点 数据 ,在 矢 
量 旋转 追踪 法 等 值 线 分 析 的 基础 上 ,经 过 滤波 、 查 找 
各 等 值 线 特性 点 .剔除 异常 点 等 处 理 ,有 序 提 了 高空 
天 气 系统 的 节点 。 

我 国 北方 处 在 中 纬度 地 区 , 主要 受 西风 带 系 统 
影响 , 且 干 时 的 发 生 往往 是 处 在 大 槽 大 疹 的 环流 转 
换 中 ,同时 由 于 西风 带 系统 的 移动 均 相对 比较 稳定 、 
有 序 , 利 于 开展 识别 实验 。 本 文 尝试 采用 计算 机 图 
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邵 建 等 :基于 图 形 学 的 致 旱 天 气 系统 自 动 识别 技术 


形 识别 技术 ,开展 致 旱 天 气 系统 的 自动 识别 ,通过 编 
程 实验 ,实现 了 高 空 槽 线 、 浓 线 、 高 低压 中 心 等 主要 
天 气 系统 的 识别 技术 ,为 天 气 预报 客观 自动 化 的 实 
现 提 供 了 一 条 新 途径 。 


1 数据 格式 说 明 


本 文 针 对 2000 一 2017 年 的 38 次 区 域 性 气象 干 
旱 过 程 ( 依 据 气 象 干旱 国家 标准 “确定 的 干旱 过 


程 ) ,利用 过 程 发 生 期 间 的 MICAPS 格 点 资料 .NECP 
逐日 再 分 析 资料 输出 格 点 资料 ,进行 系统 识别 试验 。 

MICAPS 常规 资料 包括 各 时 次 观测 实况 和 
T213 、T639 .欧洲 中 心 等 数值 预报 产品 ,其 中 风 场 资 
料 采用 第 二 类 和 第 十 一 类 MICAPS 数据 格式 ,其 他 
物理 量 采用 第 四 类 MICAPS 数据 格式 。NECP 逐日 
再 分 析 资 料 为 日 平均 基本 要 素 场 再 分 析 格 点 数据 ， 
为 “NC” 格 式 ( 表 1)。 


表 1 各 类 资料 格式 说 明 
Tab.1 Description of the data format 


资料 种 类 格 点 范围 分 辨 率 时 间 精 度 
实况 高 空 (32° ~160°E,12° ~ 80°N) 高 空 4° x4° 高 空 :08 :00 ,20 :00 
地 面 (75° ~135°E,25° ~65°N) 地 面 2° x2° 地 面 : 逐 3 bh, 以 20:00 为 主 
T213 30° ~179°E,0° ~90°N 1° x1° 每 日 20:00 起 报 , 逐 24 bh 间隔 
T639 0° ~180°E,0° ~90°N 1° x1° 每 日 20:00 起 报 , 逐 24 bh 间隔 
欧洲 中 心 30° ~179°E,0° ~90°N 2 .5x2 5 每 日 20:00 起 报 , 逐 24 h 间隔 
NECP 全 球 ,选取 30° ~179°E,0° ~90°N 2.5° x2.5° 每 日 更 新 1 次 , 逐 24 h 间隔 


注 : 由 于 2010 年 前 的 数据 保存 条 件 所 限 ,多 为 T213 、 欧 洲 中 心 和 NECP 数据 ,2010 年 后 开始 为 T639 、 欧 洲 中 心 和 NECP 数据 ;实况 数据 中 所 用 


时 刻 不 固定 ,多 以 20:00 为 主 ,无 20:00 资料 则 采用 离 20 :00 最 近 时 刻 的 数据 。 


2 天 气 系统 识别 技术 


2.1 等 值 线 曲 率 计算 

对 等 值 线条 或 系统 线条 进行 曲率 计算 ,可 得 到 
线条 的 凸凹 分 布 状况 。 一 般 的 高 度 场 中 , 槽 位 于 吓 
曲线 处 , 疹 位 于 四 曲线 处 。 而 对 于 高 空 槽 来 说 ,目的 
程度 越 强 , 则 表明 系统 中 夹带 的 冷 空 气 越 强 。 因 此 ， 
针对 高 度 场 的 曲率 计算 ,可 能 找到 平 直 气 流 ( 西北 
气流 ) 中 夹带 的 弱 冷 空气 ;而 通过 定义 曲率 强度 ,又 


可 得 到 槽 着 的 凹凸 程度, 进而 判断 槽 疹 移 速 、 变 
化 等 。 
曲率 的 计算 公式 为 : 
(XY -wy ) (1) 


[Kx 十 (y’ | 
式 中 :zy 分 别 为 xy 的 一 阶 、 二 阶 导 数 。 由 
于 曲率 大 的 量 级 平均 达 10- ,不 易 比 较 , 因 此 定义 
曲率 强度 来 进行 比较 。 
曲率 强度 Sx 定义 为 : 
Sx -hE (2) 


式 中 :6 为 曲线 中 的 最 大 曲率 ;为 该 曲线 的 平均 
曲率 。5, 为 负 时 ,曲线 向 南 ( 或 西 ) 凸 出 ;5, 为 正 时 ， 
曲线 向 北 ( 或 东 ) 凸 出 。Sx 的 绝对 值 越 大 ,表明 曲线 
的 四 凸 度 越 强 ,反之 表明 曲线 的 四 凸 度 越 弱 。 


2.2 高 空 树 脊 的 识别 

依据 模 线 定义 , 槽 点 是 某 条 等 值 线 中 曲率 最 大 
或 是 位 势 高 度 最 低 的 点 (一 般 对 西风 环流 来 讲 是 最 
南端 的 点 ) 。 分 析 覃 线 是 要 分 析 每 条 等 值 线 上 的 曲 
率 相对 最 高 点 或 位 势 高 度 相 对 最 低 点 (笔者 定义 为 
槽 点 ) , 则 可 能 出 现 某 条 等 值 线 上 存在 一 个 或 几 个 
缘 点 。 按 照 不 同 线 值 的 等 值 线 自 西向 东 分 析 ( 图 
1a) ,寻找 不 同 层 次 高 空 槽 的 槽 点 位 置 ,形成 槽 点 点 
集 ; 按 一 定 规则 连接 各 节点 则 成 槽 线 。 另 外 ,闭合 型 
等 值 线 和 开放 型 等 值 线 也 有 区 别 。 需 要 注意 的 是 不 


> 


能 让 任何 一 条 槽 线 穿 过 闭合 等 值 线 的 两 端 , 也 就 是 
Ca) 槽 线 p68 
C 1 等 高 线 
二 线 
但 、 ) a Pa 等 值 线 节点 
a 模 点 
~ \ 一 ee 
名 OD 等 高 线 走向 
Se 576 等 高 线 线 值 
Da。 ___ 着 线 
pn 等 值 线 节点 
pa es ets 
一 等 高 线 
568 f/f // 等 高 线 走向 
572 7 576 等 高 线 线 值 
576 
图 1 模 线 、 疹 线 识 别 示 意图 


Fig.1 Recognition of trough and ridge 


于 时 


说 不 能 有 一 条 槽 线 同 时 与 闭合 等 值 线 有 两 个 交点 。 

采用 类 似 槽 点 的 追踪 算法 ,对 位 势 高 度 场 和 高 
空 风 场 格 点 进行 分 析 , 得 到 高 空 状 的 脊 节点 及 状 
(图 1b)。 

2.3 系统 中 心 

系统 中 心 只 与 闭合 型 等 值 线 有 关 。 找 到 每 条 闭 
合 型 等 值 线 中 的 最 高 和 最 低 点 以 及 与 其 相对 应 的 数 
据 , 将 其 按 顺序 记录 到 数组 中 。 然 后 分 析 它 们 的 等 
值 线 值 ,找到 最 大 和 最 小 的 点 ,用 其 相对 应 数组 中 的 
数据 求 出 x 坐标 和 y 坐标 的 平均 值 , 即 为 系统 中 心 
的 坐标 。 并 且 做 好 标记 ,高 (上 暧 ) 中 心 标 记 “G(N)”， 
低 ( 冷 ) 中 心 标记 “D(L)”。 

在 目前 的 识别 过 程 中 笔者 发 现 , 当 使 用 MI- 
CAPS 格 点 资料 时 , 由 于 数值 预报 格 点 场 中 数据 分 
布 原因 ,往往 会 出 现 很 多 小 的 高 低压 系统 ”7]， 
因此 需要 改进 算法 ,找到 主要 的 天 气 系统 。 

改进 技术 的 主要 思路 如 下 : 

(1) 根据 某 一 判断 间 值 ,在 阔 值 内 出 现 的 所 有 
高 压 中 心 (或 低压 中 心 ) 计算 重心 ,其 重心 的 位 置 可 
认为 是 高 压 中 心 ( 低 压 中 心 ) 的 位 置 ,平均 强度 认为 
是 该 高 压 ( 低 压 ) 的 强度 。 

(2) 根据 某 一 判断 间 值 , 取 阔 值 内 出 现 的 所 有 
高 压 中 心 ( 或 低压 中 心 ) 中 的 最 大 值 ( 最 小 值 ) 所 在 
位 置 为 高 压 ( 低 压 ) 中 心 位 置 , 而 高 压 ( 低压 ) 强 度 即 
为 最 大 (最 小 ) 值 。 

(3) 不 考虑 几何 分 布 问 题 , 直接 取 阔 值 内 高 压 
(低压 ) 中 心平 均 位 置 为 高 压 ( 低 压 ) 中 心 位 置 ,强度 
取 阅 值 内 高 压 平均 强度 。 


3 ”关键 技术 处 理 


3.1 短 系 统 处 理 

当 分 析 识 别 出 的 系统 节点 仅 有 2 个 时 , 称 为 短 
系统 。 若 直接 连接 节点 , 则 短 系 统 为 一 条 直线 ,而 纯 
粹 直线 的 天 气 系统 几乎 是 不 存在 的 。 因 此 ,需要 根 
据 预 报 员 经 验 对 短 系统 进行 自动 化 修正 。 

按照 系统 延伸 方向 ,以 黄金 分 割 比例 增加 必要 
的 节点 ,并 用 样 条 函数 平滑 ,达到 修正 的 日 的 (图 
2)。 连 接 2 节点 构成 长 度 为 的 线段 。 若 增加 N 
点 ,过 线段 (V+1) 均 分 点 画 垂 线 , 向 某 一 方向 (根据 
正常 槽 脊 情况 , 槽 向 东 或 向 南 , 疹 向 西 或 向 北 ) 延 伸 
出 合适 长 度 ( 经 实验 ,可 取 r/4 长 度 ) ,得 到 新 增 点 。 
增加 点 数 越 多 ,系统 曲线 越 平 滑 ,但 计算 量 也 越 大 。 


一 一 ~ 槽 ( 脊 ) 线 
”| 系统 节点 
新 增 节点 
垂直 辅助 线 


图 2 短 系 统 处 理 示意 图 


Fig.2 Post-processing of short weather system 


3.2 闭合 等 值 线 长 短 轴 及 偏心 率 的 计算 
对 于 闭合 等 值 线 如 地 面 高 奈 、 高 空冷 涡 等 ,其 移 
动 发 展 与 其 形状 有 关 。 依 据 天 气动 力学 原理 ， 
对 于 闭合 天 气 系统 来 说 ,其 移动 方向 往往 与 其 长 轴 
指向 方向 重合 ,同时 与 其 重心 所 在 方向 一 致 。 因 此 ， 
计算 其 长 轴 \ 短 轴 的 指向 及 偏心 率 ,对 推测 系统 的 移 
动 有 帮助 。 

长 短 轴 的 计算 :依次 计算 闭合 等 值 线 中 每 个 节 
点 相互 间 的 距离 , 则 可 准确 地 找 出 长 轴 短 轴 。 

偏心 率 的 计算 :对 标准 椭圆 来 说 ,椭圆 两 焦点 间 
距离 和 长 轴 长 度 的 比值 称 为 偏心 率 。 但 实际 的 闭合 
曲线 往往 不 是 规则 形状 ,因此 定义 其 短 轴 与 长 轴 的 
比 为 偏心 率 。 这 也 能 反映 闭合 曲线 的 一 些 特征 。 


4 系统 识别 误差 情况 


针对 2000 一 2017 年 的 38 次 区 域 性 气象 干旱 过 
程 ,利用 改进 后 的 系统 识别 技术 进行 高 空 模 、 消 、 急 
流 轴 的 识别 检验 ( 表 2)。 按 照 系 统 的 尺度 分 别 进行 
检验 ,为 了 简便 运算 ,规定 模 线 南北 端点 距离 宇 10 
个 纬 距 ( 或 识别 节点 数 三 5 ) 的 档 为 长 波 模 , < 10 个 
纬 中 (或 识别 节点 数 <5) 的 为 短波 柳 。 

从 表 2 来 看 ,该 识别 技术 对 高 低压 中 心 的 识别 
效果 较 好 ,识别 率 均 高 于 95% ,平均 误差 都 小 于 
0.3"。 针 对 高 空 槽 线 的 识别 则 良 劳 不 齐 ,其 中 对 长 
波 槽 的 识别 精度 较 高 ,识别 准确 率 达 87% 以 上 , 平 
均 误 差 为 0.73°; 而 对 于 短波 系统 的 识别 精度 略 低 ， 
识别 率 为 74. 6% ,平均 误差 达 1. 42 个 经 纬度 。 初 
步 分 析 认 为 是 由 于 长 波 槽 系统 节点 多 ,在 系统 衔接 
和 逻辑 判断 中 具有 优势 ,更 易 判 别 ;由 于 短波 系统 节 
点 少 ,逻辑 判断 上 存在 一 定 的 难度 , 较 难 判别 。 

利用 风 场 对 识别 结果 订正 后 ,可 以 发 现 各 系统 
的 识别 准确 率 都 有 较为 明显 的 提升 、 误 差 均 有 所 减 
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表 2 主要 系统 识别 精度 检验 


Tab.2 Test of recognition precision of main weather systems 


类 型 平均 误差 /(? ) 最 小 误差 /(°) 最 大 误差 /(°) 识别 准确 率 /% 
追踪 法 风 订 正 追踪 法 风 订 正 追踪 法 风 订 正 追踪 法 风 订 正 
长 波 权 0.98 0.73 0.21 0.10 1.8 1.8 85.6 87.3 
短波 权 1,72 1.42 0.32 0.12 4.0 3.9 71.3 74.6 
低压 中 心 0. 45 0. 25 0.12 0. 10 1.47 l.2 90.3 98.5 
高 压 中 心 0.34 0.22 0.21 0.17 1.54 1.5 91.6 95.3 
少 ,尤其 高 低压 中 心 的 识别 准确 率 提高 幅度 最 为 明 流量 的 影响 [J]. 干旱 区 研究 ,2014,31 (1) :125 - 130. [ Wang 
品 ,其 中 低压 中 心 识 别 准 确 率 高 达 98. 5% ， 高 压 中 心 Zhijie, Wang Xinhui, Dan Murenjiapu, et al. Effect of climate 
se change on annual runoff volume of the Tailan River (J]. Arid Zone 
为 95. 3%; 而 长 波 模 和 短波 醒 分 别提 高 1. 7% 和 Research ,2014 ,31 (1) :125 - 130. ] 
3.3% 。 高 低压 中 心平 均 识 别 误差 减 小 到 0.3° 以 下 ， [4] 陈 伙 龙 , 王 怡 苦 , 吴 泽 旗 ,等 . 气候 变化 和 人 类 活动 对 干旱 区 内 
而 长 波 梭 和 短波 酸 的 识别 误差 缩小 到 0.73" 和 1.42。。 陆 河 径流 量 的 影响 一 一 以 新 疆 玛 纳 斯 河流 域 肯 斯 瓦特 水 文 站 
为 例 [ 圈 .干旱 区 研究 ,2015 ,32(4):692 -697. [Chen Fulong， 
S 小 结 Wang 0 Zebin ,et al. peo lia 0 and 
man activities on runoff of continental river in arid areas: Taking 
Kensiwate hydrological station in Xinjiang Manas River Basin as an 
( 1 ) 利用 追踪 算法 分 析 等 值 线 可 以 较 好 地 反应 example[ J]. Arid Zone Research ,2015 ,32(4) :692 -697. ] 
等 值 线 走 向 和 曲率 ,具有 效果 好 ,速度 快 的 特点 。 这 [5] 王 婷 , 石 育 中 , 刘 倩 ,等 . 陀 中 半 干 旱 区 不 同类 型 农户 抗旱 适应 
种 方法 在 国内 还 属于 初步 尝试 。 从 较 常用 的 统计 识 能 力 定量 评价 一 以 检 中 县 为 例 []]. 干旱 区 研究 ,2018 ,35 
别 法 .模糊 识别 法 、 基 于 特征 量 的 涡 旋 识别 等 方 nl 
Quantitative evaluation of adaptive capacity to drought resistance of 
法 ™ ww. 来 说 ? 其 算法 更 加 简 便 ? 运 算 代 人 码 编 译 更 为 different types of farmers in semiarid area in central Gansu Prov- 
快捷 。 同时 本 文 所 述 方法 不 是 对 图 像 进 行 处 理 识 ince:A case study of Yuzhong County [J]. Arid Zone Research ， 
别 , 而 是 誉 和 用 图 形 学 的 判 识 理念 和 原理 ,对 原始 2039 io a 
、 上 AAA er A EES 6] 郑 艺 , 张 丽 , 周 宇 , 等 . 1982 一 2012 年 全 球 干 旱 区 植被 变化 及 驱 
 ， 2 0 动因 子 分 析 [J]. 干旱 区 研究 ,2017,34 (1):59 - 66. [Zheng 
J 、 肖 艳 婉 可 .后 平 J \ 李 永 癌 等 所 研究 Yi,Zhang Li,Zhou Yu,et al. Vegetation change and its driving fac- 
的 图 像 识 别 技术 存在 本 质 上 的 不 同 吕 tors in global drylands during the period of 1982 ~2012[(J]. Arid 
(2) 依据 天 气 学 定义 ,利用 追踪 算法 识别 高 空 。 2 Reweeh2017 3401) 59 -时 
ee 7] 肯 大 川 , 姚 文 艺 , 焦 鹏 ,等 . 黄河 上 游 头 道 扣 站 年 径流 答 沙 系 现 
的 偏差 ;利用 风 场 进行 订正 后 ， 准确 度 明 显 提高 o 在 936.[ Ran Dachuan, Yao Wenyi, Jiao Peng, et al. Identification 
主要 系统 的 识 别 中 高 低压 中 人 心 的 识 别 玖 高 5 误差 and comprehensive diagnosis of sharp change of annual runoff vol- 
小 5 而 长 波 槽 的 识别 效果 次 之 短波 本 的 识别 效果 较 ume and silt discharge series at Toudaoguai hydrometric station in 
差 ,误差 也 较 大。 但 整体 识别 效果 基本 满足 西北 地 ee . ee ee 
区 的 气象 干旱 业务 运行 标准 。 [8] 胡 文 东 , 赵 光平 , 丁 建 军 ,等 .宁夏 短 时 强 对 流 灾害 性 天 气 预 报 
业务 工作 平台 系统 开发 策略 与 技术 措施 []]. 干旱 区 资源 与 环 
参考 文献 ( References) 和 境 ,2006,20(6):114 - 117. [Hu Wendong, Zhao Guangping, 
Ding Jianjun , et al. Develop strategies and technical measures of 
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Based on the pattern recognition technique and weather systems which could cause drought were tested 


with tracking algorithm and vector analysis methods by using MICAPS, NCEP and EC-Thin data. The results showed 


that :CD Synoptic scale troughs and ridges could be recognized with tracking algorithm ,and their curvature character- 


istics could also be obtained; ©) By using these methods ,the recognition rates of pressure center could be as high as 


90% ,and the mean errors were lower than 0.5°. However, the recognition rates of troughs varied in a range of 


74% —87% ,and the mean errors in a range of 0.7° —1.4°. With revising the wind field ,the accuracy was obviously 


improved ,and the errors were reduced to some extent ,which could basically meet the needs of predicting drought. 
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